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Abstract. Along with some fundamental concepts regarding imbibition of liquids in 
porous hosts we present an experimental, gravimetric study on the capillarity-driven 

invasion dynamics of water and of the rod-like liquid crystal octyloxycyanobiphenyl 
(80CB) in networks of pores a few nanometers across in monolithic silica glass (Vycor). 
We observe, in agreement with theoretical predictions, square-root of time invasion 
dynamics and a sticky velocity boundary condition for both liquids investigated. 
Temperaturc-dcpcHidc^nt spontaneous imbibition experiments on 80CB reveal the 
existence of a paraneinatic phase due to the molecular alignment induced by the pore 
walls even at temperatures well beyond the clearing point. The ever present velocity 
gradient in the pores is likely to further enhance this ordering phenomenon and prevent 
any layering in molecular stacks, eventually resulting in a suppression of the smectic 
in favor of the nematic phase. 



PACS numbers: 47.55.nb, 47.61.-k, 87.19.rh, 47.57.Lj, 47.56.+r 
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1. Introduction 

Liquid flow propelled by capillary forces is one of the most important transport 
mechanisms in porous environments. It is governed by a fascinating interplay of 
interfacial, viscous drag as well as gravitational forces which liquids encounter upon 
invasion into geometries with often complex topologies, such as capillary networks of 
trees (l}]3) or interconnected fractures in soils and ice. There has been signiflcant 
progress with regard to a quantitative understanding of this phenomenon, most 
prominently with regard to the transport through rocks and soils, starting with the 
seminal work of the French engineer Henry Darcy |4||5]. Almost all experimental and 
theoretical studies, published to date, deal however with porous structures where the 
characteristic length scales, meaning the typical pore diameters, are on the macroscale, 
that is much larger than the typical size of the basic building blocks (molecules) of the 
liquids |6|. 

After an introduction of the theoretical and experimental fundamentals, exemplifled 
for the case of water imbibition in mesopores, we present a temperature-dependent 
imbibition study on the liquid crystal octyloxycyanobiphenyl (80CB). The goal is to 
scrutinize which principles regarding its rheology, known from the bulk state, can be 
safely transferred to the transport across nano- or mesoscopic geometries. Finally, 
insights into the phase transition behavior of the spatially nano-conflned, rod-like liquid 
crystal will be gained. 

Given the emerging interest in micro- and nanofluidic applications [T] such a study 
on the non-equilibrium behavior of liquids in extreme spatial conflnement is not only 
of fundamental interest, but also of rather practical importance |8 10 . Moreover, the 



increasing use of mesoporous matrices as hard templates 11 for the preparation of 
well-deflned nanoscopic, soft-matter structures, such as nanotubes and nanorods 12 , 
necessitates a profound understanding of this phenomenology. 

2. Fundamentals of liquid imbibition in mesoporous solids 

2.1. Liquid flow in isotropic pore networks: Darcy 's law 

In order to introduce the fundamentals of liquid flow in mesoporous media, we start 
with a much simpler phenomenology: The flow of a liquid through a single straight 
pipe. According to the law of Hagen-Poiseuille the volume flow rate Vhp for a given 
pressure difference Ap applied along the cylindrical duct with radius r and length i is 
determined by 

I'HP^f^Ap, (1) 

Here r) denotes the dynamic viscosity of the liquid. In the following we will present a 
simple concept that allows one to extend the applicability of Eq. ([T]) to the flow through 
a complex pore network [l]. 
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Porous Vycor glass (code 7930) provided by Corning Incorporated was used over 
the course of this study. It is produced through metastable phase separation in an 
alkali-borosilicate system followed by an extraction of the alkali-rich phase. The final 
glass consists of a sponge-like network of tortuous and interconnected pores with mean 



pore radii on the nanometer scale embedded in a matrix mostly consisting of Si02 13 

Such an isotropic pore network can be characterized by a set of three quantities. 
The pore radius tq and the volume porosity 0o are probably the most intuitive ones 
among them. Two different batches of Vycor glass were applied in this study that 
differed in tq whereas they coincided in 0o = 0.32 ± 0.01. We will refer to them as V5 



[ro = (3.4 ± 0.1) nm] and VIO [tq = (4.9 ± 0.1) nm] from now on 14 . 

These matrix properties were accurately ascertained employing nitrogen sorption 
isotherms conducted at 77 K and through a subsequent analysis within a mean field 



model for capillary condensation in nano- and mesopores 15 . We obtained pore size 
distributions P{r) with standard deviations of approximately 0.1 tq (with the peak 
or most probable value of the distribution tq). In order to assess the impact of this 
uncertainty in r we turn to the quantity measured in our experiments that is the mass 
increase of the porous matrix due to the uptake of a liquid - see section |3] for details. 
Here each pore gives a signal that is proportional to its volume oc and consequently the 
macroscopic signal, which is averaged over the complete sample, should be identically 
equal to the result for an imaginary porous host with the single pore radius 



/ P(r) dr 

\ I P{r) dr ■ ^ ' 

Since we could not find any difference between tq and r' exceeding their error margins 
we conclude that the distributions can be considered to be narrow and negligible in 
these specific measurements. Hence we proceed taking into consideration only the most 
probable pore radii tq stated above. 

To account for the isotropy of the network the so-called tortuosity r must be 
introduced along with the transformation V — )• V . Pores totally aligned in flow 
direction yield r = 1, whereas isotropically distributed pores result in r = 3 since for 
randomly oriented channels only every third pore is subjected to the pressure gradient 
and hence contributes to the flow. To date many different methods have been applied 



in order to extract the tortuosity of porous Vycor glass 16 -18 . They consistently yield 
r 3.6 which, interestingly, deviates from r = 3. This is because the pores are not 
straight but rather meandering. A value r > 3 corrects the pore length for the larger 
flow path. 

Finally one is able to derive an expression that describes the flow of a liquid through 
a porous network. For a given porous matrix with cross-sectional area A and thickness 
d (along which the pressure drop Ap is applied) the volume flow rate Vd is determined 
by Darcy's law |19!| 

y^ = ^P^Ap- (3) 
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2.2. Velocity boundary conditions 



The law of Hagen-Poiseuille implies the no-slip boundary condition. This means that the 
velocity of the fluid layers directly adjacent to the restricting walls equal the velocity 
of the walls themselves. Nowadays it is indisputable that this assumption does not 
hold unreservedly. Already 60 years ago Peter Debye and Robert Cleland introduced 
both slipping and sticking fluid layers at the pore walls in order to interpret a seminal 



experiment on liquid flow across porous Vycor 19 



Up to date many factors have been found that seem to influence the boundary 
The least controversially discussed amongst them is the fluid-wall 
Shear rates beyond a critical value are supposed to induce slip, 



conditions 20 22 



interaction 23 34 



too 35 38 . In contrast, the influence of surface roughness is rather debatable 24 39 41 



Furthermore, dissolved gases 20 ,42 , the shape of the fluid molecules 28 or the add-on 



of surfactants 43 might influence the boundary conditions as well. 




Figure 1: Volume flow rate V{b) normalized with Darcy's flow rate Vb Eq. (|3j) as a 
function of the slip length b in units of the pore radius tq: x = ^. Note that the flow 
rate vanishes for b — )■ — ro -v^ x — ?■ — 1. 



The degree of slip can be quantified by the slip length b with = tq + b the 
hydrodynamic pore radius, which measures the distance from the pore center to the 
radius where the streaming velocity reaches zero. Thus, a sticking layer boundary 
condition corresponds to 6 < 0. The volume flow rate V{b) can be calculated from the 
parabolic velocity profile 

Ap 



Ar]> 



(ro + by - 



(4) 



through an integration over the channel's cross-sectional area ttTq if b > 0. For 6 < 
however the immobile layers have to be excluded. Then V{b) in units of Vd is given by 

Vi^^Ux + l)' if 6<0 

Vb \ 1 + 4x + 2x2 if 6 > ^ ^ 

in terms of normalized slip lengths x = ^. In Fig. Tj this quantity is plottet for 



— ro < 6 < ro. Because of the obvious bijectivity of V{b) any measured flow rate V 
can unambiguously be related to a unique slip length b. Part of the evaluation of the 
present investigation is based upon this principle. 
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2. 3. Dynamics of spontaneous imbibition 

From the physisist's point of view spontaneous imbibition is an impressive example for 
interfacial physics. The driving force behind the capillary rise process is the Laplace 
pressure Pl acting on the curved meniscus of a liquid in a pore or porous structure. It 
is specified by 

2 a cos 6*0 

PL = (6) 

ro 

with a the surface tension of the liquid. This implies that spontaneous imbibition can 
only occur for a wetting or partially wetting liquid. 

As the liquid rises beyond its bulk reservoir to a certain level h the hydrostatic 
pressure ph = pgh acting on the liquid column increases. The final state is derived 



from a balance between pl and ph 44 resulting in the maximum rise level given by 
Jurin's law. However, because of the tiny pore diameters considered in our study typical 
Laplace pressures are on the order of a few 100 bar and hence much higher than typical 
hydrostatic pressures. The gravitational force may therefore be neglected. 

In order to gain information on the dynamics of the imbibition process in a porous 
network we can now resort to Darcy's law. Considering Fig. [2] one easily concludes that 
the sample height d appearing in Eq. ^ has to be replaced by the time-dependent 
rise level h{t) since only the filled parts of the sample contribute to the flow dynamics. 
Moreover, at any given time t both the imbibed fluid volume V{t) and its mass m{t) 
are closely related to h(t) (utilizing the fluid's density p) via 

V{t) = '^ = <PoAh{t). (7) 

Eventually this leads to the simple differential equation 

h{t) h{t) = Ap (8) 
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solved through 



hit) 



.-^ApVt 



(9) 



the Lucas- Washburn (LW) law 45 ,46 . Substituting the Laplace pressure for Ap and 
utilizing Eq. ^ finally yields 



m(t) = p 



2rr/ 



V~t 



(10) 



c 



the mass increase m{t) of the sample due to the liquid uptake as a function of the 



time:}:. It is noteworthy that molecular dynamics simulations 47,48 corroborate on 



a microscopic scale the validity of this model for imbibition of molecular fiuids as 
derived from macroscopic, rather phenomenological considerations even down to meso- 
and nanopores. 

The imbibition coefficient C is a simple measure of the invasion dynamics. Please 
note that the above equation explicitly obeys the no-slip boundary condition. Hence, 
any deviation of a measured imbibition coefficient from the one predicted according to 
Eq. ( 10 ) may directly be interpreted in terms of a non-zero slip length b utilizing Eq. (fsl). 



3. Experimental Part 

Imbibition dynamics are studied via a recording of the samples' mass increase m(t). 
The setup is depicted in Fig. |3j For a time-dependent measurement of the force 
acting on and, hence, of the mass increase of the porous Vycor block the sample 
is installed on standard laboratory scales applying a special mounting. In order to 
perform measurements beyond room temperature a cell was constructed that allows for 
a simultaneous thermostatting of the sample itself and the liquid reservoir beneath. 

Upon start of the experiment the cell is moved upward until the sample touches the 
liquid surface. The liquid immediately starts invading the porous host what becomes 
noticeable by the gradual increase of the sample's mass (see Fig. 111). According to the 



previously introduced model of imbibition and Eq. (10) this behavior can be very well 
described by a y/t-fit (see solid line in Fig. [Z]), which yields the imbibition coefficient C. 
Eventually a level of saturation is entered, which signals that the sample is completely 
filled. 



3.1. Data interpretation using the example of mesopore- confined water 

Experiments on the capillary rise of water in V5 and VIO were conducted at three 
different temperatures T (25 °C, 40 °C, and 60 °C). The measured flow rates showed 

X Since all liquids considered in our study totally wet the high-energy silica surface of the pore walls 
(as confirmed by contact angle measurements in the group of Karin Jacobs at Saarland University, 
Saarbruecken, Germany) the contact angle is already eliminated because of cos(0°) = 1. 
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Figure 3: Exploded view of the imbibition setup (raytracing illustration). The top 
of the Vycor sample is glued to a wire and hangs freely into the cell. The cell itself 
consists of a casing and a container. Both are build out of copper and can separately 
be thermostatted. The cell itself can be moved in vertical direction. 
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Figure 4: Measurement of the mass increase of a porous Vycor block upon the liquid 



invasion during an imbibition experiment. The solid line is a A/t-fit according to Eq. ( 10 ) 
The data density is reduced by a factor of 2000. 



significant deviations from the expected value according to Eq. (10). Utilizing Eq. (|5]) 
this discrepancy could be transferred into a non-zero slip length. Consistently for both 
sample batches a negative slip length of about 5 A was found. Considering the diameter 
of a water molecule of approximately 2.5 A this result can be interpreted as follows: 
two layers of water directly adjacent to the pore walls are immobile meaning that they 
are pinned and do not take part in the flow. The residual inner compartment of the 
liquid however obeys classical hydrodynamics in accordance with former findings on the 
conserved fluidity 49 51 and capillarity 52 53 of confined water. 

The idea of the compartmentation of mesopore-confined water is corroborated by 
recent molecular dynamics studies on the glassy structure of water boundary layers 
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in Vycor and the expected existence of sticky boundary layers in Hagen-Poiseuille 



nanochannel flows for strong fluid- wall interactions 54 59 . By means of X-ray 
diffraction distortions of the hydrogen-bonded network of water near silica surfaces 
were found |60j, which might be responsible for the markedly altered liquid properties. 



Tip-surface measurements document a sudden increase in the viscosity by orders of 



magnitude in 0.5 nm proximity to hydrophilic glass surfaces 61 63 . It also extends 
former experimental results with respect to the validity of the no-slip boundary condition 
for water/silica interfaces, in which this condition was proven down to at least 10 nm 



from the surface 64 whereas slip-flow for water is only expected at hydrophobic 



surfaces 64 , 65 . For more detailed information we would like to refer the interested 



reader to two publications reporting on capillary rise dynamics of water in mesoporous 



silica 66 67 



It is interesting to note that such a compartmentisation of the pore fluid with respect 
to its dynamics has also been inferred for the linear n-alkanes confined in mesopores, 



both from imbibition experiments 68 and from neutron-spin echo experiments on the 
self-diffusion dynamics (691. 



4. Rheology and phase transition behavior of the liquid crystal 80CB in 
mesopore-confinement 

Usually phase transitions are accompanied by unique variations of the fluid properties. 
Since the capillary rise dynamics, and in particular the imbibition coefficient C, 
sensitively depend on these quantities one can easily detect such characteristic deviations 
by performing measurements at different temperatures in the T-range of interest. It is 
convenient to normalize the extracted prefactors C(T) with the value = C{T^) at 
an arbitrarily chosen temperature Tn. This procedure eliminates all quantities that do 
not depend on the temperature. Particularly, the geometry (A) as well as the internal 
topology of the substrate (ro, 0o, and r) do not play a role any longer. Only the fluid 
properties cr(T), ri{T), and p{T) enter this quantity which will henceforth be referred to 
as normalized imbibition speed Va{T) 



^(T) . r^(T,) ^^^^ 



One should not be confused by the term 'speed', since is actually dimensionless. 
Nevertheless, it quantifies the mass uptake rate of the sample at a given temperature 
and a characteristic variation in any fluid property must be mapped in this quantity. 

Here we present a thorough study on the capillary rise dynamics of the thermotropic 
liquid crystal octyloxycyanobiphenyl (80CB) in mesopores. The building blocks of this 
liquid may be considered to be rigid rods with a length of ~ 2 nm and diameter between 
Inm and 1.5 nm. Bulk 80CB undergoes a transition from the smectic A to the smectic 
N phase at T = 67 °C and beyond the clearing point = 80 °C any orientational 
ordering of the molecules disappears (isotropic phase). In particular this last transition 
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Figure 5: The experimental Couette flow conditions for measurements of the three 
Miesowicz shear viscosity coefficients of nematic liquid crystals: rji for n J^v and f? 1 1 V f , 
rj2 for f? 1 1 -y and ri ± V f , 773 for n A^v and n ± V f . 



is accompanied by the occurrence of a characteristic shear viscosity minimum, what 
might uniquely influence the invasion dynamics. In the following this behavior will be 
described in detail. 

4.1. Nemato-Hydrodynamics: Shear viscosity minimum and presmectic divergence of 
flowing nematic liquid crystals 

The momentum transport in nematic liquid crystals shows an anisotropy since it depends 
on the mutual orientations of the macroscopic molecular alignment (the director n), the 
flow velocity {v) and the velocity gradient (Vf). In 1946 Miesowicz defined three 
principal shear viscosity coefficients of nematics |70|, which can be measured in three 
different Couette flow experiments sketched in Fig. |5| Typically magnetic fields are 
applied in order to align the molecules in the nematic sample. Intuition suggests that 
the lowest resistance to the nematic flow, i.e. the lowest viscosity value, should be 772- 
Among the two remaining viscosities, rji should have the highest value. 

Indeed, when the orienting magnetic field, i.e. the director n, is parallel to the 
velocity v of the nematic flow, the lowest viscosity is recorded (see open symbols in 
Fig. [g]). Nevertheless, this relatively simple picture of the viscosity of nematic liquid 
crystals is disturbed for the compounds exhibiting the transition to the smectic A phase. 
Then, with decreasing temperature, the viscosity 172 shows a strong increase and goes 
to infinity at the N-A phase transition. The t]i and 773 viscosities are almost unaffected. 
At a temperature that is a few degrees below the temperature at which the transition 
to the smectic A phase takes place, the viscosities r]2 and r/3 interchange their roles and 
then the lowest nematic viscosity corresponds to the flow in the 773 configuration. The 
presmectic behavior of the 7/2 viscosity is due to the formation of precursors of smectic 
planes with z\\v that would be immediately destroyed by the velocity gradient, thus 
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Figure 6: Miesowicz shear viscosities ri2 (□) and r^s (A) of the hquid crystal 80CB 
compared to its free flow viscosity ?7iso ( • ) according to 71 . 



this configuration is rendered unfavorable. 

The behavior of the freely flowing compound obeys a general principle that can 
be formulated as follows: a free fluid adopts such a manner of flow, as corresponds to 
the minimum of its viscosity at given conditions |72|. Accordingly, the transition from 
the isotropic to the nematic phase manifests itself in a strong decrease of the shear 
viscosity ?7iso that is very close to ri2. Consistently, beyond the presmectic cross-over 
of ?72 and rj^ the viscosity of the freely flowing liquid crystal r^jso follows rj^, (see filled 
symbols in Fig. |6|. This result is interpreted in terms of rearrangements of the molecular 
alignment n with respect to the velocity field v, which can easily be assessed by means 
of examinations of the compound's viscosity. 



4-2. Results and interpretation 

In Fig. [rjfive measurements of 80CB invading V5 are shown along with y^-fits. Their 
coincidence is evident and hence the solution of the differential equation Eq. ^ through 
a -y/t-law is justified. This implicitly means that the right-hand side of Eq. ^ is not 
a function of the time at all, what has only tacidly been assumed so far. Especially 
the viscosity consequently is not a function of neither the time nor the time-dependent 
shear rate, which is related to the rise speed oc Shear thinning or thickening 
effects of the liquid crystal can therefore be excluded. 

The extracted normalized imbibition speeds Vn{T) are indicated by the single points 
in the upper panel of Fig. [sj Moreover, utilizing ?7iso(T)- values according to Fig. [6] along 
with the known T-dependencies of a and p one can also calculate Vn(T) (see solid hne 
in the upper panel of Fig. |8]). 

The comparison between measurement and prediction reveals a variety of 
astonishing features. First of all, both V5 and VIO reveal comparable characteristics of 
the invasion dynamics of the liquid crystal 80CB. Anyhow, only in the nematic phase the 
T-dependent behavior of the measured imbibition speeds coincides with the prediction 
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Figure 7: Specific mass uptake of V5 due to the imbibition of the hquid crystal 80CB 
as a function of the time for selected temperatures below and above the clearing point 
Tc ~ 80 °C. Solid lines correspond to v^-fits. The data density is reduced by a factor of 
2500. 
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Figure 8: (upper panel): Measured normalized imbibition speeds (for = 75 °C) 
of 80CB in V5 (A) and VIO (□), respectively, in comparison with values calculated 
on the basis of the viscosity values in the lower panel (— ). (lower panel): T-dependent 
viscosity of 80CB by way of comparison (see Fig. |6]for a detailed view). The dashed lines 
correspond to extrapolations of the calculated imbibition speeds and of the viscosity in 
the absence of the shear viscosity minimum and the presmectic divergence, respectively. 



based on the r/iso values presented in Fig. [6j In particular, the distinctive bump in 
the proximity of the clearing point Tc as the direct manifestation of the shear viscosity 
minimum in the theoretical behavior of is unambiguously absent. The imbibition 
speed rather increases monotonously with increasing temperature. 
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As mentioned before this distinctive feature of the viscosity at the clearing point 
is caused by the inset of an ahgnment of the molecules (in the nematic phase) with 
respect to the flow direction, that is n 1 1 {7. Therefore its absence can intuitively be 
interpreted in terms of an already existent alignment of the molecules beyond Tc. It 
is obvious to conclude that such an alignment is easily induced by the extreme spatial 
confinement to cylindrical pores with diameters that are not more than five times the 
length of the molecule itself. In addition, considering their rigidity and, in particular, 
the parabolic flow profile established in the pore, it is hard to think of any alternative 
to the tendency of molecular alignment parallel to the pore axis and, consequently, to 
the flow velocity v. From this point of view the absence of the shear viscosity minimum 
and, consistently, of the nematic to isotropic phase transition is not surprising at all 
but rather consequent. In this context it is more suitable to label the phase beyond the 
clearing point not isotropic but paranematic (P). This term preeminently discloses the 
solely confinement-induced alignment of the liquid crystal. 

It was demonstrated experimentally 73 -77 , in agreement with expectations from 



theory 78 -80 , that there is no 'true' I-N transition for liquid crystals confined in 
geometries spatially restricted in at least one direction to a few nanometers. The 
anchoring at the confining walls, quantified by a surface field, imposes a partially 
orientational, that is, a partially nematic alignment of the confined liquid crystal, even 
at temperatures T far above the clearing point Tc. The symmetry breaking does not 
occur spontaneously, as characteristic of a genuine phase transition, but is enforced over 
relevant distances by the interaction with the walls. 

The assumption of such a paranemtic phase is corroborated by recent birefringence 
measurements of rod-like liquid crystals (7CB and 80CB) confined to an array 
of parallel, nontortuous channels of 10 nm mean diameter and 300 pm length in a 



monolithic silica membrane 81,82 . These measurements elucidate that the surface 
anchoring fields render the bulk discontinuous I-N transition to a continuous P-N 
transition. The transition temperature Tc is found to be changed only marginally, 
due to a balance of its molecular alignment induced upward and its quenched disorder 
(attributable to wall irregularities) induced downward shift. This agrees with the 



observations of liquid crystals imbibed in tortuous pore networks 73 , 74 



Interestingly, due to the complete absence of the shear viscosity minimum in the 
results of the T-dependent series of imbibition measurements shown in Fig. [8] a definition 
of a P-N transition temperature in confinement is not possible at all. This is confirmed 
by a simple extrapolation of the ri2 viscosity to higher T, which saliently reproduces 
the measured imbibition speeds beyond the bulk clearing point. This can only be 
interpreted in terms of an extremely high degree of orientation already existent in the 
paranematic phase; at least higher than suggested by the birefringence measurements 



81 . This difference is most probably due to the basically differing detection method 
of the P-N transition: normally birefringence or calorimetry (DSC) measurements 
performed with the confined static liquid are applied for this purpose. But, the viscosity 
measurements presented here refer to the liquid's dynamics in mesopore confinement. 
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As already mentioned before, the additional emerging flow velocity and in particular 
the velocity gradient seemingly enhance the paranematic orientational alignment in the 
r72-configuration significantly. 

The second remarkable feature of the results presented in Fig. [8] is the absence of 
the presmectic divergence of the viscosity of the freely flowing liquid crystal. This would 
result in a dramatic drop of the imbibition speed due to the inset of smectic layering. 
However, the measured values do not indicate such an effect. Its absence rather suggests 
a suppression of the A phase in favor of the N phase. This is elucidated by a simple 
extrapolation of the viscosity to lower T as indicated in Fig. [8} which preeminently 
reproduces the measured imbibition speeds. 

What are the reasons for this discovery? First of all, the confinement to a cylindrical 
pore (rather than to a film geometry like in the Couette flow depicted in Fig. [5]) renders 
the ?73 viscosity as unfavorable as the rji viscosity. This again clarifies the high stability 
of the nematic ry2-configuration in confinement. A cross-over behavior as ascertained for 
the freely flowing bulk liquid can hence be excluded. Yet, the presmectic divergence of 
ri2 is caused by the destruction of precursors of smectic planes for n\\v. From this point 
of view the suppression of the A phase is not surprising at all but a mere consequence 
of the overall stabilization of the ?72-configuration in the mesopore confinement. 

Even the static mesopore-confined liquid crystal shows such modified mesophase 



behavior 75 . For example the heat capacity anomaly typical of the second-order 



N-A transition in rod-like liquid crystals immersed in aerogels is absent or greatly 
broadened 83 , 84 . Nuclear magnetic resonance (NMR) measurements revealed the lack 



of pretransitional smectic layering due to the rough surface of the confining walls 85 



Furthermore, a systematic study of the influence of the degree of confinement indicates 
that the N-A transition becomes progressively suppressed with decreasing pore radius 



whereas the stability range of the nematic phase is increased 82 , 86 



Finally we'd like to comment on results from an examination on the absolut 
invasion dynamics of the liquid crystal as opposed to the analysis of normalized 
imbibition coefficients just presented. This investigation was performed according to 
the procedure introduced in section 3J^ referring to mesopore-confined water. For the 
liquid crystal 80CB we determined slip lengths of 6 = (—1.11 ± 0.23) nm for V5 and 
b = (—1.54 ± 0.31) nm for VIO, respectively. Referring to the molecular shape stated 
earlier we may hence, again, conclude that one layer of flat lying molecules is pinned to 
the pore walls and does therefore not contribute to the flow dynamics. Apparent velocity 
slippage at the walls, as was expected to be associated with the molecular alignments 
in the channels 87 , could not be detected in our study. 



5. Conclusions 



After a short introduction of liquid imbibition in porous hosts we have presented a 
gravimetric study on spontaneous imbibition of water and the liquid crystal 80CB. We 
inferred from our studies an interfacial boundary layer adjacent to the pore walls with 
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a defined thickness whose dynamics are mainly determined by the interaction between 
hquid and substrate. This manifests itself in terms of a negative velocity slip length for 
both liquids investigated. The flow properties of the pore-condensed molecules in the 
pore center are, however, remarkable robust, that is bulk-like or at least similar to the 
bulk phase (in the case of the liquid crystal). 

The T-dependent imbibition measurements of the liquid crystal 80CB revealed 
that confinement plays a similar role as an external magnetic field for a spin system: 
the strong first-order I-N transition is replaced by a weak first-order or continuous 
paranematic to nematic transition, depending on the strength of the surface orientational 



field 88 . Based on detailed knowledge of the static (equilibrated) liquid's behavior in 
the mesopores as deduced from previously conducted birefringence experiments, we were 
able to procure complementary results with respect to its dynamic (non-equilibrium) 
behavior. The additional emerging flow velocity and in particular the velocity 
gradient enhance the paranematic orientational alignment significantly rendering the 
P-N transition even broader than known from the equilibrium state. Due to the high 
stabilization of the 772-configuration in the N phase the A phase is suppressed and the 
stability range of the nematic phase is increased. 

The finding of two distinct species (with regard to flow dynamics) is reminiscent 
of the partitioning of pore condensates found in vapour sorption isotherms (film- 
condensed versus capillary-condensed state) [89], of the two species often reported 
in studies on self-diffusions dynamics (slowed-down or increased self-diffusion versus 



bulk dynamics) [69 , 90 92 , of the effects of pore condensates on the deformation of 
the matrix (expansion versus contraction of the pore walls upon change from film- 
condensed to capillary-condensed state) [93]. Last but not least, it reminds of the 
distinct crystallization behaviour of pore condensates (an amorphous boundary layer 
versus a crystallized fraction of molecules in the pore center) [89]. In principle, this 
partitioning, but also the highly stabilized nematic phase of a confined liquid crystal, 
can be traced to the strong influence of the substrate potential on the liquid or solid 



layers, respectively, right adjacent to the pore walls 94 96 



From a more general perspective, our experimental studies demonstrate that 
an understanding of the non-equilibrium behavior of confined soft-matter system 
necessitates a thorough understanding of their equilibrium behavior, which quite often 



differs markedly from the macroscopic phenomenology 97 98 
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